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We present the first study of the intrinsic electrical properties of WS 2 transistors fabricated with two 
different dielectric environments WS 2 on Si0 2 and WS 2 on h-BN/Si0 2 , respectively. A comparative analysis 
of the electrical characteristics of multiple transistors fabricated from natural and synthetic WS 2 with 
various thicknesses from single- up to four-layers and over a wide temperature range from 300 K down to 
4.2 K shows that disorder intrinsic to WS 2 is currently the limiting factor of the electrical properties of this 
material. These results shed light on the role played by extrinsic factors such as charge traps in the oxide 
dielectric thought to be the cause for the commonly observed small values of charge carrier mobility in 
transition metal dichalcogenides. 



The emerging class of atomically thin semiconducting materials formed by transition metal dichalogenides 
(TMDCs) is showing a plethora of complementary properties to those of graphene that are of interest to 
fundamental and applied research. These materials are uniquely suited to study the superconducting phase 
transition in the extreme two-dimensional limit inherent to atomically thin systems 1 " 4 . At the same time TMDCs 
have a band gap which is essential for transistor applications and which could enable a new class of atomically thin 
photo-transistors. For example WS 2 has a direct band gap of 2 eV in single layer form 5 " 8 and has already shown 
great promise as a flexible transistor with field effect mobilities comparable to the best liquid crystals and on/off 
ratio of the current exceeding 10 6 9 . Understanding the limiting factors of the electrical properties of TMDCs is an 
open quest and a stepping stone for accessing novel physics in these systems. 

The typical values of charge carrier mobility measured in thin WS 2 flakes are always much lower than those 
measured in bulk material 2,4 . This behaviour has been interpreted as due to defect states in the Si0 2 substrate 
leading to the localization of charge carriers in TMDCs and a small charge carrier mobility 10 . To probe the 
intrinsic electrical properties of TMDCs it would be necessary to measure electrical transport in either suspended 
structures or in transistors fabricated on clean substrates with fewer impurities than typically present in Si0 2 . An 
ideal choice for such a substrate is hexagonal boron nitride 11 , which is a preferred substrate for high quality 
graphene transistors since it has a very low concentration of charge scattering impurities and is atomically flat 12 . 
To date such a study has not yet been conducted and the consequent lack of knowledge is limiting the potential 
impact of TMDCs on fundamental and applied research. Furthermore most of the studies conducted so far have 
been limited to just MoS 2 , while other TMDCs such as WS 2 have not yet received much attention, whereas they 
might be better suited than MoS 2 for a given application. 

Here we present the first study of the electrical properties in WS 2 transistors fabricated on different dielectrics 
(i.e. Si0 2 and h-BN/Si0 2 ) and using synthetic as well as natural WS 2 . The comparative analysis of the electrical 
characteristics of these transistors studied in the temperature range from 300 K down to 4.2 K shows that in all 
cases electrical transport takes place via hopping conduction through localized states 13,14 . At low temperature 
(T < 20 K) we observe peaks of the conductance as a function of back-gate voltage and source-drain bias due to 
inelastic tunnelling in the impurity states with sub -gap energy. These results show that intrinsic disorder rather 
than extrinsic factors such as defect states in the oxide dielectric is limiting the electrical properties of WS 2 and 
more generally TMDCs. 
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Results 

Thin flakes of WS 2 were obtained by mechanical exfoliation of flakes 
from synthetic crystals onto p-doped Si/Si0 2 substrate that serves as 
a back gate (for natural WS 2 see supporting information). Thin flakes 
are first identified with the aid of optical microscopy and their thick- 
ness is subsequently determined by atomic force microscopy (AFM) 
and Raman spectroscopy. The fabrication of WS 2 transistors on h- 
BN and subsequent encapsulation in h-BN is carried out using the 
dry transfer method first developed for graphene 12 . This consists of 
exfoliating WS 2 onto a substrate coated by water soluble polymer and 
PMMA. After dissolving in water the soluble polymer, the free WS 2 / 
PMMA bilayer is aligned onto previously exfoliated h-BN (—20 nm 
thick) on p-doped Si/Si0 2 . The substrate is then heated up to melt the 
PMMA and secure contact between WS 2 and h-BN and the PMMA is 
subsequently removed in acetone. Electrical contacts to WS 2 are 
fabricated using standard electron beam lithography, thermal evap- 
oration and lift-off of Cr/Au (5/70 nm). 

Figure la shows an AFM measurement of a thin WS 2 flake with a 
fold in the upper left corner highlighted by a dashed line. A statistical 
study of the height measured in areas which include the step edge at 
WS 2 /Si02 (region A) and the folded corner (region B) shows a 
comparable step height of —1.6 nm in A and —1.3 nm in B, see 
Figure lb. Since the thickness of a monolayer WS 2 flake is 
—0.65 nm 815 we conclude that this flake is a bilayer. A comparative 
plot of the Raman spectra (see methods) for WS 2 with different layer 
numbers shows marked differences depending on the specific thick- 
ness of the flake, see Figure lc. More specifically it is known that the 
peak with low Raman shift (—350 cm" 1 ) is a convolution of two 



Lorentzians (Figure Id) whose positions change as a function of 
the layer number 16 . One Lorentzian is due to the second order lon- 
gitudinal acoustic phonon mode (2LA(M)) corresponding to collect- 
ive oscillations of the atoms in the plane, and this gives a Raman peak 
at 352.7 cm" 1 in single layer WS 2 . The second Lorentzian is given by 
the in-plane optical phonon mode (E^fT)) representing the in- 
plane counter oscillations of W and S atoms in the lattice. Finally 
the out-of-plane optical phonon mode (A 1<? (r)) representing the 
out-of-plane oscillations of W and S atoms gives a Raman peak at 
416.6 cm" 1 . A plot of the relative wavenumber shift (Av that is the 
difference between the 2LA(M) and A 1<? (r)Raman peaks) for a large 
number of flakes with various thicknesses shows that Av changes in a 
discrete way according to the number of layers which have been 
independtly measured with AFM, see Fig. le. Finally, upon increas- 
ing the number of WS 2 layers the position of the 2LA(M) and (r) 
peaks redshift monotonously, whereas the A^(r) peak blue shifts as 
previously shown 16 , see Figure If. 

Having established a reliable procedure to identify the layer num- 
ber of WS 2 flakes we now turn to investigate the electrical transport 
properties of this material. The source -drain current vs. bias voltage 
characteristics (I-V) of WS 2 transistor devices are always highly non- 
linear and upon performing current-bias annealing, a linear I-V 
around zero voltage bias is attained (see Figure 2a and supplementary 
information). Owing to the difference in work function between WS 2 
and Cr, a Schottky barrier of about 100 meV has to be expected at 
this interface when no-gate voltage is applied. The observed bias- 
annealing changes in the I-V and the large values of voltage bias at 
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Figure 1 | (a) shows an AFM measurement of a bilayer WS 2 > the scale bar corresponds to 500 nm. The dashed areas labeled by A and B enclose the 
step edge at the Si0 2 -bilayer WS 2 and the fold in the WS 2 flake respectively. The corresponding histograms of the measured heights in A and B are shown 
in (b). (c) shows the evolution of the shape and position of the Raman peaks of WS 2 as the number of layers is increased from single layer to bulk, (d) is a 
plot of the Raman spectra for bulk WS 2 and a fit to three Lorentzians corresponding to the 2LA(M), E^fT) and A 1(? (T), see main text, (e) shows the 
measured wavenumber shift (Av) between 2LA(M) and A ig plotted for 30 flakes with different layer number, (f) summarizes the measured Raman shift 
for 2LA(M), E^CO an< ^ A lg (T) as a function of layer number. 
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Figure 2 | (a) shows plots of I-V for a 4 L WS 2 device on a 20 nm thick h-BN crystal before and after sweeoubg the voltage bias to large values (see 
supplementary information). The inset shows the circuit used for decreasing the contact resistance, (b-d) are plots of the gate dependence of the 
conductivity for a monolayer WS 2 flake on Si0 2 substrate (b), for a four layer WS 2 flake on Si0 2 substrate (c) and for a 4 layer WS 2 flake on h-BN 
substrate (d). All the sweeps in (b-d) were made at the same rate of 100 V/hr. The inset in (d) shows a micrograph picture of a WS 2 transistor with scale 
bar of 5 um. 



which these non-linearity occur suggest a different origin for this 
phenomenon, that is the possible presence of an oxide barrier at 
the WS 2 /Cr interface which can be electrically broken upon applying 
a large voltage bias as shown in Figure 2a. In the following we only 
consider the analysis of electrical transport measurements in devices 
after bias -annealing. 

Figure 2b-d show the room temperature field effect transistor 
(FET) transfer characteristics, that is the gate voltage (V^) depend- 
ence of the conductivity (cr), for monolayer WS 2 on Si0 2 (Figure 2b), 
four-layer WS 2 sample on a Si0 2 (Figure 2c) and four-layer WS 2 
sample on a h-BN/Si0 2 (Figure 2d). In all cases we observe that 
the conductivity has a large on- off ratio typical of semiconducting 
materials, with a finite threshold voltage. However we find that 
the field effect mobility (fi) is always larger in WS 2 on h-BN than 
in WS 2 on Si0 2 (0.23 cm'V-V 1 for lL-WS 2 /Si0 2 , 17 cm'V-V 1 for 
4L-WS 2 /Si0 2 and -80 cm'V-V 1 for 4L-WS 2 /h-BN/Si0 2 in 
Figure 2(b-d)). A large hysteresis is also present in <j(V g ) for WS 2 
on Si0 2 but is fully suppressed when WS 2 is on h-BN/Si0 2 . Similar 
hysteresis in I-V have also been reported in graphene and is com- 
monly attributed to dopants present in the Si0 2 dielectric 1718 . 

Discussion 

For all the measured devices we find that the temperature depend- 
ence of <r( V g ) shows a pronounced suppression of the value of a upon 
lowering the temperature as expected for a semiconducting material, 
see Figure 3a. In these devices we apply a large enough value of gate 
voltage such that the charge carriers are directly injected from the 
metal contacts into the conduction band of WS 2 . In this limit the 
relevant energy scale dominating the temperature dependence of 
the zero -bias resistance is the difference between the Fermi energy 
and the conduction band edge of the n-doped semiconductor (i.e. 
WS 2 ) 20 . A plot of a as a function of T" 1 at V g = 60.5 V reveals that 
from 260 K down to 100 K the conduction takes place by thermally 
activated charge carriers, i.e. a(T) = cr 0 exp( — 3s/2k B T) with 3s the 



activation energy and k B the Boltzman constant. The values of 3s 
estimated from a fit of o(T) for 50 V < V g < 60 V are in the range 
0.109 eV < 3s <0.113 eV and change linearly with V g , see inset in 
Figure 3b. These values of 3s are compatible with the voltage bias 
range over which non-linear I-V are measured (see blue curve in 
Figure 2a) suggesting that Ss is the energy from the Fermi level to 
the conduction band edge (E c ), i.e. 3 s = E c — E F which is also much 
larger than the Schottky barrier height (—100 meV). 

The smooth dependence of 3 s on V g demonstrates that for sub -gap 
energies the Fermi level can be continuously tuned by means of a gate 
voltage throughout the defect induced states. To estimate the density 
of defect states we consider the equivalent gate capacitance of these 
WS 2 transistors that is the series of the gate oxide capacitance (C ox ) 
and defect states capacitance (Q), i.e. — dEp/dVg = 1.5 x 10~ 4 e = 

C C 

e ox 2 . Knowing that the oxide capacitance per unit area is 

C ox + Q 

C ox = ^«1.2x 10" 4 Fm" 2 we find Q = 0.8Fm" 2 = q 2 D(E), 
a 

where q is the unit of charge and D(E) is the density of defect states 
which we estimate to be 3.12 X 10 37 J _1 m" 2 . 

The dominant role of disorder induced states with sub -gap ener- 
gies becomes fully apparent when considering a fit of the low tem- 
perature o{T) in logharitmic scale in terms of T~ p with p critical 
exponent, see Figure 3c. This study reveals that p = 1/3 gives the 
best fit stemming for non-interacting Mott variable range hop- 
ping 13 ' 1419 a = G 0 exp( — T 0 /T) V3 where T 0 is the hopping parameter 
and is related to the density of localised states existing within the 
forbidden gap and the electron wavefunction size £ by the following 
13 8 

relation T 0 = . The extracted values for the hopping para- 

k B D c £ 

meter T 0 at each different gate voltage are plotted in figure 3(d) along 
with the conductance at T = 4.2 K. This comparative plot shows a 
clear correlation between the hopping parameter and the conduc- 
tance whereby peaks in conductance correspond to very low values of 
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Figure 3 | (a) shows a representative plot of a{ V g ) for a 4L-WS 2 in the gate range 50 V to 65 V. (b) Typical temperature dependence of the conductivity 
plotted in terms of the activation energy relation (this curve is taken at V g = 60.5 V). The inset is a plot of the extracted activation energy for different gate 
voltages, each point corresponds to an average over 0.2 V gate voltage, (c) shows the same data as in (b) but plotted in terms of 2D Mott variable 
range hopping relation, (d) The conductivity at 4.2 K plotted alongside the hopping parameter T 0 . A strong correlation between the two is observed: i.e. 
peaks in conductance correspond to low values of T 0 . 

T 0 . Furthermore T 0 is found to fluctuate from —100 K to —4000 K in 
a small gate range (from V g = 50.5 V to 52 V, corresponding to an 
energy window of just 0.25 meV). Consequently the estimated local- 
ization radius in WS 2 increases from 1.8 nm to 17 nm. These obser- 
vations indicate that the sub -gap impurities states have peaks of 
narrow energy band-widths dominating electrical transport for 
sub-gap energies. 

Another prominent feature evident in the temperature depend- 
ence of (r(V g ) is the emergence of peaks for T < 100 K with decreas- 
ing amplitude for T < 20 K, see Figure 3a. At the same time the 
differential conductance as a function of source-drain bias and gate 
voltage at T = 4.2 K (Figure 4a) shows that these peaks shift their 
position as a function of voltage bias. These observations suggest that 
charge transport at sub -gap energies occurs through inhomogeneous 
charge puddles and localized states in WS 2 . Since we observe a sim- 
ilar a(V g ) behaviour in a variety of samples independently of (1) the 
WS 2 flakes aspect ratio, (2) the WS 2 layer number and (3) the dielec- 
tric environment (WS 2 /BN/Si02, see supplementary information) 
we conclude that the localized states dominating electrical transport 
in WS 2 at sub -gap energies are intrinsic to the WS 2 and not extrinsic 
such as defect states in the dielectric. 

To estimate the localization radius (£) we consider electrical trans- 
port measurements of a representative 4L-WS 2 in which the peaks of 
(j{V g ) are spaced by an average gate voltage (V^) ~ 1.13 V corres- 
ponding to 0.17 meV, see bottom graph in Figure 4. In this device the 
peaks of a at fixed V g as a function of source-drain bias ( V) are spaced 
by an (V) ~ 11 mV, which for a channel length of 350 nm corre- 
sponds to a threshold electric field E T = 3.14 X 10 4 V/m. This value 
of E T together with the observed average peak separation of 
0.17 meV gives a localization region of diameter 2£ = 5.4 nm which 
is consistent with the extracted value of the localization radius £ from 
the analysis conducted on the temperature dependence of o(V g ). 




Figure 4 | The top colour map shows the measured differential 
conductance plotted against gate voltage and source drain at T = 4.2 K 
for the same representative 4L-WS 2 on h-BN discussed in Figure 3. The 
bottom plot is a graph of the conductivity at T = 4.2 K while the inset is a 
graph of five differential conductance curves plotted from V g = 61.52 V up 
to 62 V highlighted by the dashed lines A and B respectively and in steps of 
96 mV of gate voltage. 
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Finally we note that Coulomb blockade cannot account for the 
observed peaks of o{V g ). Indeed, if we assume a charging energy in 
our devices of Ec ~ 40-50 meV estimated directly from the stability 
diagram shown in Figure 4, we extract a diameter d ~ e 2 /4e 0 8 r Ec ~ 
20-40 nm for the confining regions (e = 1.610" 19 C, £ 0 = 8.8510" 12 F/ 
m and e r = (s vac + s BN )/2 = 2.5 with the dielectric constant for 
vacuum and BN s vac = 1 and s BN = 4). Given the dimensions of 
the conductive WS 2 channel, our devices would consist of 100-1000 
charging regions (i.e. (length X width)/d = (350 nm X 1500 nm)/ 
d). The stability diagram of such an array of charging islands would 
consist of many overlapping Coulomb diamonds which are not 
observed in our measurements. An indication of the underlying 
physical process originating these peaks of a(V g ) is given by the 
temperature dependence of o{V g ) presented in Figure 3a: we alwasy 
observe that the amplitude of the peaks decreases upon lowering the 
temperature. This behaviour has been previously reported in other 
semiconducting systems 21,22 and it is a fingerprint of inelastic tun- 
nelling which in WS 2 occurs through the sub-gap impurity states. 

In summary we have presented the first systematic study of the 
intrinsic electrical properties of thin WS 2 flakes. By comparing the I- 
Vg of transistors fabricated using two different dielectric environ- 
ments (i.e. (1) WS 2 on Si0 2 and (2) WS 2 on h-BN/Si0 2 ) we find that 
hopping through localized states dominate electrical transport over a 
wide temperature range (T < 100 K). This intrinsic disorder has a 
finite density of states at sub -gap energies which contribute with 
inelastic tunnelling to electrical transport. These results demonstrate 
the dominant role played by intrinsic disorder over extrinsic factors 
such as defect states in the oxide dielectric as a limiting factor of the 
electrical properties of WS 2 . 

Methods 

Materials. Synthetic WS 2 was purchased from Lowerfriction.com. 

Measurement techniques. The Raman spectra where measured with a Renishaw 
spectrometer using an excitation laser with a wavelength of 532 nm, focused to a spot 
size of 1.5 /mi diameter and 1 mW incident power. These measurements were 
performed in air and at room temperature. 

Electrical measurements. The electrical transport measurements were performed in 
constant voltage configuration with excitation voltage smaller than k B T, with k 5 
Boltzmann constant. The differential conductance was measured using the lock-in 
technique. 
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